Barkhausen Noise in a Relaxor Ferroelectric 
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Barkhausen noise, including both periodic and aperiodic components, is found in and near the 
relaxor regime of a familiar relaxor ferroelectric, PbMg!/3Nb 2 /303, driven by a periodic electric 
field. The temperature dependences of both the amplitude and spectral form show that the size of 
the coherent dipole moment changes shrink as the relaxor regime is entered, contrary to expectations 
based on some simple models. 

PACS numbers: 77.80.Dj, 77.80.Bh, 77.84.Dy, 72.70.+m 



Although relaxor ferroelectrics form locally ferroelec- 
tric polar nanodomains, underlying chemical disorder 
somehow prevents the formation of long-range ferroelec- 
tric or antiferroelectric order Q. Instead, there is a 
crossover, characterized by faster-than-Arrhenius tem- 
perature dependence of relaxation rates (approximately 
of the Vogel-Fulcher form, e.g. Q), to a glassy "re- 
laxor" regime. The key ingredients of a successful model 
may include static random polar fields (e.g. g]), random 
anisotropy (e.g. [Q), random interactions among the po- 
lar nanodomains (e.g. 10, || ^, and interactions with 
soft TO phonons ml pi M and slow charge-transfer modes 
§• 

Nonequilibrium noise and mesoscopic noise provide 
model-sensitive probes of glassy freezing O, unlike 
macroscopic linear response and the corresponding equi- 
librium fluctuations, which show rather gen eric proper- 
ties. The Barkhausen effect (see, e.g., Q) arises be- 
cause the polarization changes unevenly as an applied 
field changes, with the random part of the response only 
statistically similar between nominally identical samples. 
For two materials with similar average response func- 
tions, the Barkhausen noise will be larger for the material 
whose polarization makes larger coherent changes. 

Although simple pictures of interacting polar nan- 
odomains would predict that the dynamically coherent 
units grow as the frozen regime is entered, a previous 
attempt to find Barkhausen noise in the relaxor regime 
found none, setting an upper limit on the cooperative 
changes of the dipole moment Large discrete steps 
in polar order have been found in the best-studied re- 
laxor, PbMg 1/3 Nb 2 /303 (PMN), but only after the ap- 
plied dc electric field was large enough to drive the ma- 
terial into an ordinary ferroelectric regime with macro- 
scopic domains j3j. Here we report Barkhausen noise in 
the nonferroelectric (paraelectric and relaxor) regimes of 
PMN, where conventional Barkhausen noise from macro- 
scopic domains is absent. 

The simplest Barkhausen model consists of a set of 
polarization steps at fixed fields distributed with uniform 



probability over the field sweep range. For periodically 
driven fields, the noise consists of a periodically repeated 
random walk, so S(f), the Fourier power spectral density 
of the voltage fluctuations, will consist only of harmonics 
of the driving frequency fp , and its spectral envelope will 
decay as l// 2 . If the step pattern changes between cycles 
(e.g., due to thermal jitter in the step times) then there 
will be aperiodic noise too. For fast-relaxing domains 
which remain in quasi-cquilibrium as the field changes, 
the polarization steps follow the smooth change in Boltz- 
mann factor vs. time (on a coarse time scale), causing the 
l// 2 envelope to cut off rapidly (csch 2 (7r///c)) above a 
characteristic frequency fc 
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(Here, Ep is the maximum amplitude of the sinusoidal 
driving field, p is the typical change in electric dipole 
moment in a single step, and kp is Boltzmann's con- 
stant.) For domains with large barriers to switching, the 
switches occur as single abrupt steps. However, ther- 
mal jitter in the timing of otherwise reproducible steps 
also reduces the periodic component above fc, with a 
(numerically computed) envelope very similar to that for 
the fast switchers. The missing harmonic power should 
appear as an aperiodic component of roughly Lorentzian 
shape and corner frequency comparable to fc- 

When this simple independent-step picture applies, the 
voltage variance (the integral of S(f)) will grow linearly 
with Ep (so long as f c > fp) P§, being 
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Here C is the geometrical capacitance, and we assume 
the dimensionless dielectric constant £> 1. When there 
is a range of p's, p in Eq. ^| is replaced by (p 2 ) / (p). 

Generically, if the non-Arrhenius slowing down of re- 
laxation rates arises because interacting nanodomains 
form clusters which reorient coherently but with barriers 
which increase with cluster size (e.g. as in the random- 
field Ising model jl4)), then the dipole moments of the 
coherent clusters (even if the constituents are randomly 
aligned) would grow as T was reduced, causing ((SV) 2 ) 
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FIG. 1: e'(T) and e"(T) at 50 Hz, measured at 13 V/cm. 
The inset shows the nonlinear out-of-phase response, most 
apparent at high-T where the linear out-of-phase response is 
tiny. 



and fc to rise. If the activation barrier heights for some 
clusters become too large for them to respond at or near 
fn for a given Ed , then the typical p would saturate near 
the moment of the largest responding clusters. The size 
of the largest responding cluster would shrink slightly as 
T is reduced further, but only as T to a power less than 
1. We found neither continued growth nor approximate 
saturation of the noise magnitude and frequency range, 
but instead abrupt shrinkage. 

Our sample consisted of two pieces, each about 
0.75 mm thick, from a single crystal of PMN, grown by 
the Chokhralsky method at the Rostov-on-Don Institute 
of Physics. The sputtered gold contact pads on the main 
(001) faces had about 1 mm 2 surface area each. 

Standard ac susceptibility measurements at 50 Hz 
showed the usual PMN relaxor behavior, as shown in 
Fig. The response was close enough to linear over the 
field range employed in the noise experiments for the lin- 
ear susceptibility to suffice for approximate calculations. 
Previous experiments Jl5| have shown that the apparent 
Vogel-Fulcher freezing temperature of PMN (w 220 K) 
determined by ac response is essentially unchanged over 
the field range used in these experiments, although it 
drops sharply at higher fields. 

For the noise measurements, the two halves of the sam- 
ple were incorporated in a balanced-bridge circuit, as 
shown in Fig. |[ This setup allowed us to apply an ac bias 
up to 10 volts on each half of the sample without over- 
loading the low-noise amplifiers, despite the generation of 
harmonics by the sample nonlinearity, slow capacitance 
changes due to systematic aging effects and imper- 
fect common-mode rejection. However, any residual sys- 
tematic differences in the nonlinearity of the two parts 
would generate harmonics. An experiment in which one 
sample arm was replaced with a polystyrene capacitor in- 
dicated significant systematic harmonics out to the fifth 
harmonic. We avoided using harmonics lower than the 
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FIG. 2: The schematic diagram of the Barkhausen measure- 
ment circuit. A trimming capacitor in the balance circuit 
has been omitted to avoid clutter. SI and S2 are the sample 
halves. ±G is the gain of the driving op-amps. Data were 
taken on the 14-bit converter in sets of 30,000 points with a 
10 kHz sampling rate, giving 1/3 Hz frequency bins in the 
Fourier spectrum. Usually 16 such spectra were averaged. 
The anti-alias filter corner frequency was 5 kHz, although it 
should have been 4 kHz. 
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FIG. 3: S(f) with E D = 130 V/cm at f D = 50 Hz at 244 K 
(spectrum 2) is shown above a background (spectrum 1) con- 
sisting of equilibrium noise from the dissipative sample plus 
amplifier noise, mainly from current noise from the Stanford 
550 preamp. 60 Hz and its harmonics have been edited out. 

tenth in the data analysis, since systematic nonlinearities 
are distinct from Barkhausen noise. The ac source has 
very low levels of noise and distortion, and the analog-to- 
digital sampling is synchronized well with the ac drive. 

Fig. H shows a typical S(f) with, roughly as expected, a 
mixture of periodic and broad aperiodic noise. The wings 
around fr> (and to a lesser extent around low harmonics) 
were not anticipated, but these can arise from small drifts 
in the bridge cancellation as the sample ages. The size of 
the polarization events, calculated below from S(f), was 
far too small for them to be discerned individually. 

Since the spectral form is close to the independent-step 
expectation, we can extract a typical step size p, although 
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FIG. 5: The magnitudes of periodic and aperiodic spectral 
components (above the relevant backgrounds) in the range 
512 to 1024 Hz, weighted by e'(T), are shown vs. T at two 
different Ed- 



FIG. 4: Noise power (above the relevant backgrounds) in 
the octave from 512 to 1024 Hz as a function of Ed for (a) 
aperiodic and (b) periodic components. 



we doubt that the polarization changes literally consist of 
simple discrete steps. As shown in Fig. fj, both periodic 
and aperiodic components of the noise grow nonlinearly 
with Ed over the range explored. Thus the p calculated 
at some Ed gives only a typical dipole step size under the 
particular drive conditions. The form of the nonlinearity 
suggests a broad distribution of p's, with many small p's 
showing up in the noise only as Ed is increased enough 
to make fc > fo for that p. The form of the aperiodic 
spectrum indicates a similar distribution. 

Although we are not yet confident in extracting abso- 
lute step sizes from this new technique (mainly due to the 
difficulty of removing all artifacts from the lower harmon- 
ics, which contain the most power), we can use Eq. | to 
calculate a lower-bound estimate for a typical p by inte- 
grating the broad aperiodic spectrum only. (We assume 
this p reflects about how big the cooperative changes are 
regardless of whether they are simple discrete steps.) At 
T = 250 K, just above the relaxor regime, at the largest 
Ed used, the result is about p = 2 x 10~ 22 C cm, which 
apparently would grow at higher Ed- Under the same 
conditions, the envelope of the harmonic noise has a l// 2 
tail extending from about 1 kHz to 4 kHz, the top of our 
measured range, indicating that there are some discrete 
dipole switches as large as 2 x 10 -21 C cm. An order- 
of-magnitude estimate of the moment of a polar nan- 
odomain at 250 K, based on the range of the static po- 
lar correlations [[L7| and the saturation polarization [[I8| , 
would be about 3 x 10" 23 C cm. Thus it seems that even 
in the paraelectric regime there are some dynamically co- 
herent units larger than single polar nanodomains. 



Despite the uncertainties in the absolute calibration, 
we can obtain a good indication of how typical p's depend 
on T, a key issue for relaxor models. For independent 
steps, the quantity e 1 (SV) 2 is proportional to the ac- 
tual p's under particular experimental conditions. Fig. ^ 
shows how some periodic and aperiodic components of 
e' (SV) 2 depend on T at fixed E D - We pick spectral com- 
ponents in an octave around / = 15/d, because in this 
range the signal-to-background is still good, we expect 
few artifacts in the harmonics, and the T-dependence of 
p shows up in S(f) both via the overall magnitude and 
via the T-dependence of fc- (Other frequency compo- 
nents show the same sort of behavior, just less clearly.) 

As seen in Fig. [|, a s the relaxor regime is approached 
from high T, the spectral components of e' (SV) 2 ei- 
ther rise slowly, as one would guess from the growth 
of medium-range static correlations jlTjl , or show little 
change. However, below 235 K to 250 K, both periodic 
and aperiodic components drop sharply, in sharp con- 
trast to the continued growth of static correlations. The 
factor e'(T) accounts for less than half of this drop. 

The width of the envelope of the harmonic components 
is also T-dependent, being largest (for Ed = 110 V/cm) 
at 250 K. In contrast to the long l// 2 tail at 250 K, 
below 250 K the harmonics fell sharply into the an- 
harmonic background, e.g. becoming undetectable above 
about 900 Hz at 210 K, confirming that the largest dipole 
moment steps occur near 250 K. 

In summary, our noise results indicate abrupt shrink- 
age of the typical net dipole moment of the dynamically 
coherent changes (at fixed driving field) as the relaxor 
regime is entered. We do not want to burden this exper- 
imental result with much theory, but some initial com- 
ments should be helpful. 

We know of no proposed model in which the coher- 
ent dipole moment differences between the equilibrium 
states at different fields would shrink sharply in the re- 
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laxor regime (i.e., relaxors are not believed to be frus- 
trated antiferroelectrics). Rather, the shrinkage of the 
steps found on a fixed frequency scale should be related to 
the Vogel-Fulcher kinetics HQ], i.e., to the rapid growth 
of barrier heights as the relaxor regime is approached. 
Although our results do not dictate a model for the re- 
laxor transition in PMN, we can rule out some models. 

A model of the dynamics with a fixed set of relaxation 
modes with a fixed broad unimodal distribution of acti- 
vation barriers (presumably monotonically increasing in 
p) iQ , would predict that the characteristic p is a weakly 
increasing function of T at all T, and hence is inconsistent 
with our results. For simple growth of independent polar 
nanodomains, the dipoles grow along with the anisotropy 
barriers, so the typical dipole moment for which the field 
(plus thermal activation) allows switches at rates of order 
fc would continue to grow. For formation of coherently 
flipping clusters, the typical moment also grows with the 
number of nanodomains (probably as a square root). Al- 
though the moment to barrier ratio would be lower for 
larger clusters, the moment of the largest clusters with 
rapid enough driven kinetics would approximately satu- 
rate. Unless there were a strange distribution of clus- 
ter sizes, so that as the typical size grew, the mean of 
those remaining below some fixed cutoff would shrink, 
such models would not reproduce our results. 

Thus the Barkhausen results strongly suggest that the 
rapid growth of barriers is not just the result of nan- 



odomain growth or of collective nanodomain clusters 
forming. Instead, it seems that each nanodomain's con- 
tribution to the barriers grows much more rapidly than 
its contribution to the dipole moment, as in a recent 
phenomenological model [Q, suggesting that less-polar 
degrees of freedom are freezing along with the polar nan- 
odomain orientations. Another indication of some such 
effects is that the low-temperature linear heat capacity 
specifically associated with the relaxor [^0| gives a di- 
mensionless entropy of about 10 3 per nanodomain. 

Models of relaxor freezing invoking strong coupling of 
nanodomain polarization to multiple slow local degrees 
of freedom (particularly charge transfers) and to soft TO 
phonon modes |§| have been proposed to account for clus- 
ters with multiple metastable polarization states having 
a roughly uniform distribution in solid angle ||. Such 
models would allow complex aging |l6[ ] even in small 
clusters of nanodomains. Whether or not such mod- 
els prove successful, the Barkhausen data provide a new 
qualitative constraint on models for the relaxor transi- 
tion, which has thus far provided a target for somewhat 
under-constrained theory. 
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